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A detailed thermo-economic model considering different technological alternatives for 
thermochemical production of Synthetic Natural Gas (SNG) from lignocellulosic biomass is 
presented. First, candidate technology for processes based on biomass gasification and 
subsequent methanation is discussed and assembled in a general superstructure. Both 
energetic and economic models for biomass drying with air or steam, thermal pretreat¬ 
ment by torrefaction or pyrolysis, indirectly and directly heated gasification, methane 
synthesis and carbon dioxide removal by physical absorption, pressure swing adsorption 
and polymeric membranes are then developed. Performance computations for the 
different process steps and some exemplary technology scenarios of integrated plants are 
carried out, and overall energy and exergy efficiencies in the range of 69-76% and 63-69%, 
respectively, are assessed. For these scenarios, the production cost of SNG including the 
investment depreciation is estimated to 76-107 € MWh _1 SNG for a plant capacity of 
20 MW thi biomass> whereas 59-97 € MWh _1 SNG might be reached at scales of 150 MW th) biomass 
and above. Based on this work, a future thermo-economic optimisation will allow for 
determining the most promising options for the polygeneration of fuel, power and heat. 

© 2009 Elsevier Ltd. All rights reserved. 


1. Introduction 

1.1. Background and goal 

Biomass is a widely available, renewable energy source whose 
use is generally accepted as a contribution to decrease 
greenhouse gas emissions by substituting fossil fuels. In 
Switzerland, the total forest area is increasing and the 
sustainable use of wood as an energy source could at least be 
doubled in the medium term, reaching a value of 
12.1 TWhyear -1 [1]. Its share in the total energy consumption 
is currently of 2.6% and it is almost exclusively used as a fuel 
for heating purposes. In the same time, the emissions of the 


transport sector account for 42% (16.5 Mtons) of the total fossil 
C0 2 emissions. Despite the efforts to mitigate these emissions 
according to the Kyoto protocol, they are expected to increase 
[2] and still few solutions to enhance the use of renewable 
fuels are emerging. In this context, processes that convert 
lignocellulosic biomass into Synthetic Natural Gas (SNG), 
which is equal in quality to fossil natural gas, are currently 
under investigation [3-6]. Such a product could easily be 
injected into the gas grid and thus profit from the existing 
distribution network for transport applications. Along with an 
increasing market share of gas engines in the transport sector, 
fossil fuels could thus be partially substituted by a renewable 
fuel that is neutral in greenhouse gas emissions. Furthermore, 
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Nomenclature 

kc 2 /C! 

proportionality constant c G 2 h 4 /cch 4 , - 



I 

length, m 

Abbreviations 

YYl 

mass flow, kg/s 

CFB 

circulating fluidised bed 

n 

stoichiometric coefficient, mol 

FICFB 

fast internally circulating fluidised bed 

P 

pressure, bar 

PSA 

pressure swing adsorption 

Q 

heat, kW 

SNG 

synthetic natural gas 

R 

ideal gas constant, kJ/(K mol) 

TSA 

temperature swing adsorption 

SN 

stoichiometric number, - 

vol 

volume 

T 

temperature, K 

wt 

weight 

trl 

relative duration of PSA adsorption, - 

Greek letters 

Ah 0 lower heating value, kj/kg 

dh° standard heat of reaction, kj/mol 

Ah uap heat of vaporisation, kj/kg 

tr2 

U 

U 

V 

relative duration of PSA recycling, - 
velocity, m/s 

overall transfer coefficient, variable 
volume flow, m 3 /s 

Ak° 

exergy value, MJ/kg 

Subscripts 

£ 

energy efficiency, % 

BM 

bare module 

V 

exergy efficiency, % 

cc 

carbon conversion 

£ 

mass efficiency, % 

eg 

cold gas 


humidity, kg„ 20/kgtot 

d 

drying 

0 

relative humidity, % 

el 

electric 

P 

density, kg/m 3 

eq 

equilibrium 

Roman 

A 

letters 

area, m 2 

GR 

9 

Im 

m 

sat 

th 

tot 

vap 

grass roots 
gasification 
log-mean 
methanation 

saturation 

thermal 

total 

vaporisation 

b 

C 

cost exponent, - 
cost, € or €/MWh 

c 

molar fraction, % 

d 

diameter, m 2 

E 

e 

mechanical/electrical power, kW 
specific mechanical/electrical work, kj/kg 

9 

molar Gibbs free energy, kj/mol 

Superscripts 

G 

gas mass flux, kg/(s-m 2 ) 

+ 

Material or energy stream entering the system 

h 

mass enthalpy, kj/kg 

- 

Material or energy stream leaving the system 

h 

height, m 

0 

Standard conditions 

K 

equilibrium constant, variable 




as carbon dioxide has to be separated from the methane as 
a byproduct, a negative C0 2 balance would be obtained if the 
C0 2 is sequestrated. 

Recent studies indicate that thermochemical SNG produc¬ 
tion by means of biomass gasification is technically feasible 
and energetic efficiencies ranging from 58% to 70% have been 
assessed [3-7]. Economic analysis has further shown that 
the overall production costs are expected to amount to 
28.1 €MWh _1 for a plant capacity of 100MW th based on the 
lower heating value of wood [3]. Similar research on biomass 
gasification and subsequent Fischer-Tropsch synthesis claims 
that liquid fuels could be produced at conversion efficiencies of 
33-50% and at a cost of 57.0 € MWh -1 for a first commercial 
plant having a nominal input capacity of 367 MW th [8]. The 
exergetic efficiency of such facilities might reach 46% [9]. With 
total energy efficiencies of around 50% [10,11], the production 
of ethanol from lignocellulosic biomass performs similar to 
Fischer-Tropsch products, but the co-production of SNG 
instead of power from the residual lignin can increase the 
biomass to fuel conversion to 70% [11]. Most of these conclu¬ 
sions are drawn from investigations of the existing technology 
for the subsequent process steps and considering specific 


scenarios. If modelled at all, the chemical transformations 
of the process streams have mostly been calculated assuming 
thermodynamic equilibrium, whereas it is known that this 
is not accurate in case of gasification [12]. Furthermore, only 
limited energy integrations have been performed, cogenera¬ 
tion possibilities have not been studied in detail and the 
process economics have been analysed with capacity-based 
correlations disregarding the specific process conditions. 

The present work is based on the investigation of the 
energetic performance considering a specific process design 
without proper process integration [5]. Its purpose is to 
develop detailed thermodynamic and economic models of the 
technological alternatives for thermochemical SNG produc¬ 
tion using woody biomass as raw material. In particular, it is 
aimed to identify the optimal topology and operating condi¬ 
tions of the process with regard to efficiency and cost. To do 
so, an approach based on the definition of a process super¬ 
structure and using advanced process integration techniques 
as well as multi-objective optimisation has been applied 
[13,14]. This paper deals with the identification of candidate 
technology and it’s modelling in order to prepare a future 
thermo-economic process optimisation. 
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1.2. Methodology 

The present paper follows a previously developed method¬ 
ology addressing the optimal thermo-economic process 
design for thermochemical fuel production from biomass [15]. 
According to this approach, a preliminary analysis of the 
process specifications and the required steps is carried out to 
identify suitable technology for the process. A process 
superstructure including all possible options is then defined. 
For each unit, thermodynamic and economic models are 
developed. The thermodynamic model consists of two parts, 
i.e. energy-flow models representing the chemical and 
physical transformations of the chemical species from raw 
materials into products and an energy-integration model that 
represents the heat recovery in the system. In the energy- 
flow models, mass- and energy balances as well as physical 
properties and compositions of all process streams are 
computed using the commercial flowsheet calculation soft¬ 
ware Belsim-Vali [16]. The energy requirements associated 
with these transformations define the hot and cold streams 
that are integrated in the energy-integration model using the 
heat cascade concept. In this model, utility streams whose 
flow rates are determined to satisfy the process units 
requirements are considered together with the process 
streams. By solving a linear programming problem, the 
optimal process layout and its corresponding utility system is 
determined with regard to operating cost [17]. The data from 
these thermodynamic models are further used to estimate 
the capital costs of the installation. Considering specific 
process conditions like temperatures, pressures, volume 
flows and gas compositions, available data from existing 
plants are used for a preliminary sizing of the major process 
equipment. The costs of the single components are then 
estimated from correlations [18,19] and the total grass roots 
cost of the installation is determined [19]. 

The major advantage of constructing the process model 
with separate energy-flow and energy-integration models is 
that it is not necessary to define a priori a specific process 
topology and heat exchanger network like it is the case in 
conventional process simulation based on a scenario 
approach. With this approach, it is possible to consider a 
maximum number of potential solutions without setting any 
preferences already at an early stage of the design. Contrary, 
by coupling the model to an optimisation software both 
advantageous operating conditions of the process and its 
optimal layout are simultaneously determined and sub- 
optimal solutions due to unnecessary restrictions on the 
topology are prevented. Therefore, the methodology is well 
suited for process design and has been applied successfully 
to fuel cell system [13] and power plant design [14]. 


2. Process description 

2.1. Process block flow 

The general layout of the thermochemical process converting 
lignocellulosic biomass to SNG consist of three major steps, 
namely gasification (endothermic), methane synthesis 
(exothermic) and gas separation. Representing the feedstock 


by a typical wood composition and using the carbon atom as 
reference, the overall conversion is exothermic and can be 
expressed as: 

Ah r =-10.5 kJxmoD+d 

CHi 35O0 63 T 0.3475 H 2 0 0 0.51125 CH 4 (1) 

+0.48875 C0 2 

The considered raw material is wood chips whose proper¬ 
ties are given in Table 1. Due to its high moisture content, 
a drying stage prior to gasification is necessary. Furthermore, 
the gas produced through gasification needs to be cleaned 
from impurities to prevent catalyst damage during methane 
synthesis. As biomass contains too much oxygen to be 
completely reformed into methane, C0 2 is by-produced and 
must be removed before or after methanation in order to meet 
the quality requirements of natural gas. In Switzerland, 
a Wobbe Index of 13.3-15.7 kWh Nm' 3 and a methane content 
of at least 96%vol are required for unlimited injection into the 
national high grid that operates at around 50 bar [20]. 

2.2. Technology for the process 

2.2.1. Wood drying 

The high moisture content of wood at the gasifier inlet 
severely decreases its performance [12]. This is mainly 
because high-temperature heat from above the process pinch 
is consumed for water desorption and evaporation in the 
gasifier, which appears as an important exergy loss [22]. The 
presence of steam further tends to shift the equilibrium of the 
gas phase towards higher H 2 and C0 2 fractions and decreases 
the CH 4 content in the gas, which has also been observed 
experimentally by for example Gil et al. [23]. A drying stage 
before gasification is thus essential. Its level is subject to 
optimisation since steam is also used as gasifying agent. 

Steam and air drying are reported to be the most common 
technologies applied in sawdust or wood chips treatment [24- 
26]. The main difference between these processes is that the 
use of steam allows to efficiently recover the consumed heat 
by condensing the produced additional steam at its boiling 
temperature, while this energy is normally lost in air drying. 
However, operating temperatures are usually higher in steam 
dryers, and air drying might be advantageous if heat is avail¬ 
able at lower temperatures [25,26]. The performance of both 
processes depends on the integration with the rest of the 
process and its heat recovery possibilities. The optimal choice 
is thus determined by energy-integration aspects, and models 
for both steam and air drying have been developed. 

2.2.2. Gasification 

Gasification of wood is an endothermal process where solid 
macromolecules are broken into mainly hydrogen, carbon 


Table 1 - Proximate and ultimate analysis of wood fuel. 


Proximate analysis 


Ultimate analysis 

Ah 0 

18.6 MJ/kgdry 

C 

51.09%wt 

Ak 0a 

20.9 MJ/kgdry 

H 

5.75%wt 

Humidity (4>) 

50.0%wt 

O 

42.97%wt 

Ash content 

0.6%wt 

N 

0.19%wt 


a Chemical exergy has been calculated using the method of [21]. 










1590 


BIOMASS AND BIOENERGY 33 (2OO9) 1587-1604 


monoxide, carbon dioxide, hydrocarbons, tars and ash. 
According to the requirement for a final gas product with 
high calorific value, adequate gasification technology should 
produce a nitrogen-free gas with high methane content. Air as 
gasifying agent is therefore not suitable and only steam and 
oxygen can be used for this purpose. Considering the equilib¬ 
rium equations for the gas phase [5], high methane fractions 
are expected for gasification at low temperature and high 
pressure, having the further advantage that endothermicity of 
gasification decreases. Entrained flow gasification technology 
operating at high temperature is thus not adequate for the 
targeted application. The specified plant capacity further 
restricts the gasifier choice. Due to geometric considerations, 
fixed bed reactors are limited to about 10 MW th [3,27] and thus 
hardly compatible with plants at industrial scale. 

In this study, two types of gasifiers have been investigated, 
i.e. indirectly heated, steam blown FICFB-gasification [28] 
operating at around 850 °C and atmospheric pressure and 
directly heated, steam-oxygen blown, pressurised CFB-gasifi- 
cation [29] operating at around 800 °C. 

2.2.3. Oxygen production for gasification 

Oxygen required for gasification is conventionally produced 
by cryogenic air separation or pressure swing adsorption 
(PSA). According to Kirschner [30], on-site production with 
these technologies gets competitive at capacities of about 
1000 Nm 3 h -1 (0.35 kgs -1 ) and 50Nm 3 h -1 (0.02 kgs -1 ), 

respectively. If smaller flow rates are required, oxygen can be 
purchased and delivered to the plant from an external 
supplier. Apart from these established technologies, high- 
temperature air separation by ceramic ion transfer 
membranes may be promising in the future since they can 
benefit from a tight integration into the plant [31]. If tempo¬ 
rarily cheap electricity is available, electrolysis is further an 
interesting option since the co-produced hydrogen can be 
injected into the methane synthesis, where it is bound to the 
abundant carbon from biomass and increases the SNG-yield. 
Even if not operated as base load, this would allow for peak 
shaving electricity generated from intermittent sources (like 
wind power) and store it as green fuel in the gas grid [32]. 

2.2.4. Thermal pretreatment before gasification 

In addition to drying, the biomass feed can optionally be 
processed in a second, thermochemical torrefaction or 
pyrolysis pretreatment step. In the literature, these technol¬ 
ogies are often discussed in the context of improving the solid 
fuel’s thermal and mechanical properties like heating value or 
grindability [33,34], or the direct production of bio-oil that is 


further refined to liquid fuel [35]. However, thermochemical 
processes based on gasification can benefit in general from 
such a pretreatment [22,33]. In addition to completely dry the 
feedstock for gasification, these technologies are charac¬ 
terised by an onsetting endothermal decomposition at low 
temperature, which decreases the energy demand at high- 
temperature in the subsequent gasification. This directly 
results in a higher cold gas efficiency if the required energy for 
torrefaction or pyrolysis is provided from excess heat below 
the pinch. The product gas of a directly heated gasifier is thus 
less oxidised, and the fuel consumption of an indirectly 
heated gasifier is reduced. 

As demonstrated on pilot scale by Henriksen et al. [36], one 
alternative is to directly close-couple a pyrolysis screw with the 
gasification stage and thus feed both the gaseous and solid 
products into the gasifier. Another one is to only feed the solid 
product to the gasifier and burn the by-produced low calorific 
gas. Although Prins et al. [33] concluded that this option 
decreases the performance of directly heated fluidised bed 
gasification, it is promising in indirectly heated gasification 
since the volatiles can be used as fuel. 

2.2.5. Gas cleaning 

During gasification, tars are formed and traces of nitrogen, 
sulphur, chlorines and metals contained in the wood are 
reformed and transferred into the product gas. In order to 
prevent catalyst poisoning, the gas must be rid of these 
substances prior to methanation. Conventional cold gas 
cleaning includes a baghouse or sand filter to remove solid 
particles and partially tars, a scrubber for removal of 
ammonia, metals and residual tars as well as guard beds for 
scavenging hydrogen sulfide. Typical temperatures for these 
stages are 150-180 °C at the filter inlet, 40 °C at the scrubber 
outlet and around 350 °C in the guard beds [37,38]. Alterna¬ 
tively, hot gas cleaning by particle removal with candle filters 
or electrostatic precipitators, thermal or catalytic cracking of 
the tars and high-temperature adsorption of other contami¬ 
nants could be applied. This would allow for a compact 
process design based on pressurised gasification and metha¬ 
nation without intermediate gas cooling and compression. A 
general overview on these advanced cleaning technologies for 
biomass gasification is given in for example [39], extensive 
details on catalytic processes are reviewed by [40] and recent 
technology developments can be found in [41-43]. 

2.2.6. Methane synthesis 

Methane synthesis is a refining process to increase the calo¬ 
rific value of a gas containing high carbon monoxide and 


Table 2 - Gasification and methanation reactions. 

Interaction 

Name 

Reaction 


Jh° 

Solid-gas 

Hydrogenating gasification 

C(s) + 2H 2 ^CH 4 

( 3 ) 

-75 kj/mol 


Boudouard equilibrium 

C(s) + C0 2 ^2C0 

( 4 ) 

173 kj/mol 

Gas-gas 

Methane synthesis 

CO + 3H 2 ^ CH 4 + H 2 0 

( 5 ) 

-206 kj/mol 


Ethene reforming 

C 2 H 4 + 2H 2 0 ^ 2CO + 4H 2 

(6) 

210 kj/mol 


Water-gas shift equilibrium 

co + h 2 o^±co 2 + h 2 

( 7 ) 

-41 kj/mol 
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hydrogen fractions. The principal conversion is described by 
the methanation reaction (5) in Table 2. At typical operating 
temperatures of 300-400 °C and preferably under pressure, 
higher hydrocarbons are broken down to CO and H 2 and form 
additional CH 4 , as shown for ethene (6). Depending on the 
initial gas composition, carbon dioxide is further reformed or 
produced through the water-gas shift equilibrium (7). 

The stoichiometric coefficients of Equations (5)-(7) allow to 
determine the amount of hydrogen that is needed to 
completely reform CO, C0 2 and C 2 H 4 into methane. It is thus 
convenient to define the stoichiometric number SN of the 
incoming gas stream in order to characterise the achievable 
methane yield [44], which becomes for the considered 
reactions: 


SN ch 4 = 


Ch 2 

3c C o + 4 cco 2 + 2cc 2 h 4 


( 2 ) 


To obtain a highly pure methane stream, this ratio must be 
close to unity. If it is below, the feed gas lacks hydrogen and 
the product gas will contain a non-negligible amount of 
carbon dioxide, which is the case for gases originating from 
wood gasification. 

Methane synthesis of a H 2 /C0/C0 2 mixture is highly 
exothermic and its reactor design is critical with regard to 
temperature control. Common installations use product gas 
recycle loops or multiple intercooled reactors with prior steam 
addition [44,45]. Alternatively, an internally cooled fluidised 
bed reactor for isothermal once-through methanation has 
been developed and successfully operated at pilot scale for gas 
produced by coal gasification [46]. Ongoing research has 
further proven the suitability of the latter to gases produced by 
wood gasification [38]. 


2.2.7. Carbon dioxide removal 

Removal of C0 2 from natural gas is a standard operation in gas 
refining applications. The choice of the optimal technology 
depends on specific process conditions like the amount of gas 
treated, the partial pressure of carbon dioxide in the feed and 
the required purity of the produced gas [47]. For gases con¬ 
taining high C0 2 fractions of more than 40% as it is the case in 
this application, physical absorption [48], pressure swing 
adsorption (PSA) [49] and membrane processes [50] can be 
identified as the most suitable options 1 . Optimal operating 
pressures of these processes depend considerably on the 
pressure of the raw gas. In general, physical absorption and 
membrane processes typically adapt to the gas grid pressure 
of 50 bar, while the maximum pressure in a PSA cycle is of 5- 
6 bar. Contrary to chemical absorption considered by [6], all 
this processes are almost neutral in thermal energy. 


2.3. Process superstructure 


the utility system including hot and cold utilities as well as 
cogeneration possibilities by a steam Rankine cycle is not 
detailed. In addition to the mandatory operations outlined in 
Section 2.1, the superstructure includes the optional thermal 
pretreatment and also the possibility to adjust the stoichi¬ 
ometry of the producer gas prior to methanation. Disregarded 
by [3-5,38], this is a popular option for coal-derived producer 
gas that is recently also considered for biomass [6,51]. 


3. Process modelling 

3.1. Energy-flow models 

As detailed in the applied methodology, the energy-flow 
models aim at defining the chemical and physical trans¬ 
formations occurring in the units of the superstructure. They 
determine the heat transfer requirements in terms of heat 
load and temperatures to be satisfied by the heat recovery and 
utility system, and further provide the necessary data for the 
equipment rating. In this regard, the process design method 
implies a trade-off between the details of the unit models and 
the degrees of freedom left to the process design. It requires 
models that are simple and robust enough to allow for eval¬ 
uating different process configurations, and precise enough to 
accurately represent the performance of the process units. 
The proposed models realise this compromise without 
requiring detailed simulations of the heat- and mass-transfer 
dynamics and the chemical reaction kinetics. This approach is 
demonstrated in the following sections, and general 
assumptions and parameters are shown in Table 3. 

3.1.1. Drying 

Air drying processes are governed by the transfer of water in 
the solid and differences in vapour pressure at its surface and 
the surrounding air. The equilibrium vapour pressure at the 
solid’s surface depends on temperature and is given by sorp¬ 
tion isotherms. To quantify this equilibrium, an empiric 
equation has been fitted to data from [53]: 

0 wood = 2.865-10 -VJr + (2.307-10- 1 

- 1.273-10 _3 (T - 273 ))(p air - 2.519-1CT 1 <?4 
+ (2.199-1CT 1 + 8.630-10 _4 (T — 273))^ 3 ir (8) 

In this equation, <£^ 0 ^ terms the wood humidity in 
kg H2 o-kg~l, (p a i r the relative humidity of air in % and T the air 
temperature in K. Limited residence times in dryers prevent 
however that thermodynamic equilibrium is reached, and 
mass transfer needs to be modelled in order to estimate the 
outlet moisture content. In analogy to the energy transfer 
equation in heat exchanger calculations, the amount of 
transferred moisture can be written as: 


From the preceding considerations, a general process super¬ 
structure is identified and depicted in Fig. 1. The internal 
subsystem configuration like recycling and the integration of 


1 Although chemical absorption with amines is a widely used 
technology for acid gas removal, it is not considered here since it 
is better suited for feeds at lower partial pressure of C0 2 . 


^H 2 Q, uap 

hlair 


DpApi m 


(9) 


where U p terms an overall mass-transfer coefficient and Ap Im 
the countercurrent log-mean difference of partial pressures, 
i.e.: 


Apim 


Api - A p 2 
In (Ap a /Ap 2 ) 


( 10 ) 
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Fig. 1 - Process superstructure including main process streams without recycling loops. Dashed lines assemble investigated 
alternatives for different process sections and dotted lines indicate optional units. 


In this equation, Api terms the difference of partial pressures 
at the air outlet and Ap 2 at the air inlet. U p has been fitted to 
data of a direct rotary drum dryer for wood [24] and a value of 
11.16-10 -3 bar -1 has been determined. 

Contrary to air dryers, the gas phase in steam dryers is 
homogeneous and always above the moisture’s boiling 
temperature. The drying process is thus controlled by heat 
transfer. As for air drying, an overall transfer coefficient U T 
based on the amount of energy needed to evaporate the 
moisture has been introduced: 


Ah uap m H2 o, ua p 

team 


— U T ATi m 


( 11 ) 


U T has been fitted to data for a pilot-scale steam dryer [25] and 
a value of 1117 J kg -1 K -1 has been determined. 

The performance of both dryers is evaluated in terms of 
drying efficiency e d and the specific mechanical energy 
consumption e + d . The former is defined as the ratio between 
the energy used for moisture evaporation and the total 
thermal energy supplied to the dryer (Eq. (12)) and the latter as 
the mechanical work consumed per unit mass of evaporated 
moisture (Eq. (13)): 


£d 


Ah ua pm H2 o, uap 

Qd 


( 12 ) 


e 



(13) 


Fig. 2 compares the computed performance of the drying 
processes as a function of the temperature of the heating fluid 
at the dryer inlet. In case of steam drying, e d is slightly 
decreasing with temperature and situated at around 60%. The 
performance of air drying is expected to increase with 
temperature and reaches efficiencies above 70% at high 
temperatures. For both systems, the specific mechanical 
energy consumption is decreasing with temperature and is in 
the order of 200 kj kg H20 -i Mp . 


3.1.2. Gasification 

The thermochemical decomposition of solid feedstock is 
a complex process that is governed by chemical kinetics and 
heat- and mass-transfer dynamics. As discussed by for 
example Pierucci and Ranzi [54], accurate models to predict 
the performance of gasifiers require to handle the description 
of a large number of species and reactions, and thus to 
numerically solve a large system of differential equations. 
Such models, however, are not appropriate for process design 
purposes. In flowsheeting applications, the product compo¬ 
sition and yield is therefore often simply fixed or computed 
through an equilibrium approach. Although sometimes 
assumed (e.g. [12]), thermodynamic equilibrium is yet not 
a valid assumption for fixed and fluid bed gasification that are 
typically operating below 1000 °C, and various corrections 
have been proposed. Many authors, among recent studies [55- 
57], therefore fix the carbon conversion and/or the fraction of 
methane and higher hydrocarbons, or correct the equilibrium 
constant with a multiplication factor [58,59]. Introduced by 
[60] and discussed for its application to biomass gasification 



Fig. 2 - Drying efficiency and specific mechanical energy 
consumption for steam and air drying. 
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by for instance [22,61] and in more detail by [62-64], a ther¬ 
modynamically more significant approach is to correct the 
equilibrium temperature by introducing a temperature 
difference AT eq to equilibrium, i.e.: 


K p =K p (T g + AT eq ) (14) 

where K p is the theoretical equilibrium constant and K p the 
apparent one corresponding tom the experimentally observed 
composition at the gasification temperature T g . They are 
computed with: 


/EM“( T )-E n i9?( T A 

K P (T)=exp 2-- 

V / 


Up? 


(15) 


(16) 


where n are the stoichiometric coefficients of the reactants i 
and products j, p their partial pressure, g the Gibbs free energy 
and R the ideal gas constant. Assuming that wood is converted 
into the gaseous components H 2 , CO, C0 2 , CH 4 , C 2 H 4 , H 2 0, N 2 
and residual solid carbon C(s), four model equations are 
needed in addition to the atomic balances in order to deter¬ 
mine the product yield and composition. Among the different 
possibilities, the hydrogenating gasification (3), Boudouard (4), 


steam ethene reforming (6) and water-gas shift (7) equilibria 
of Table 2 could be chosen as an independent set. 

According to the heat of reaction of these equilibria, higher 
gasification temperatures would favour the formation of H 2 
and CO, while C0 2 , CH 4 , C 2 H 4 and - with the typical wood 
composition of Table 1 - C(s) are favoured at lower tempera¬ 
tures. This behaviour is generally confirmed for H 2 , CH 4 , C 2 H 4 
and C(s) by some extensive studies on fluidised bed gasifica¬ 
tion [23,55,65]. In indirectly heated gasification, this trend has 
also been roughly confirmed for CO and C0 2 [65], whereas 
directly heated gasification with oxygen generated more C0 2 
and less CO with increasing temperature [23]. These differ¬ 
ences can partly be explained by the increasing oxygen 
content with temperature due to the direct heating, but must 
also be attributed to higher reaction rates that favour the 
conversion towards C0 2 . 

A large amount of different operating conditions has been 
investigated in the literature, yet surprisingly few information 
is available on the influence of pressure. To our knowledge, no 
systematic dataset has been published, and the few 
measurements at moderate and only slightly varying pressure 
do not allow to confirm the expectation from Eq. (16) that 
namely the CH 4 -fraction increases with increasing pressure. 
Indeed, it can be argued that the pressure tends to accelerate 
the reaction rates, thus improves the conversion towards 
equilibrium and counterbalances the positive equilibrium 
effect of pressure on methane yield. We therefore adopt 


Table 3 - General assumptions and key parameters for the energy-flow models. 1 

Section 

Specification 

Value 

Air drying 

Wood outlet temperature 

Tair.out 


Dryer pressure drop 

100 mbar 

Steam drying 

Wood outlet temperature 

TH 2 0,sat 


Dryer pressure drop 

50 mbar 


Heat loss (based on energy for steam heating) 

18% 

Gasification 

Pressure drop in fluidised bed reactors 

150 mbar 


Excess pressure of injected steam 

12 bar 

Gasification, indirectly heated 

Reactor heat loss (based on transferred heat) 

10% 


N 2 -content in cold gas 

0.5%vol 

Gasification, directly heated 

Reactor heat loss (based on Ah° wood ) 

1% 

Gas cleaning 

Filter inlet temperature 

150 °C 


Pressure drop 

100 mbar 


Biodiesel consumption in scrubber 

4.7 ml-m“ 3 g as 

Methane synthesis 

Reactor pressure drop 

150 mbar 


Energy efficiency of electrolysis 

85% 

Physical absorption 

C0 2 -solubility in Selexol 

0.18 moll -1 bar -1 


Relative solubility C0 2 /CH 4 

17.1 


Regeneration pressure 

1 bar 


Pressure drop through column 

neglected 

Pressure swing absorption 

Adsorption pressure 

5.5 bar 


Purging pressure 

0.2 bar 


C0 2 (CH 4 )-slip fraction during adsorption (purging) 

0.2%vol 

Polymeric membranes 

Material: cellulose acetate (polysulfone) for bulk-C0 2 (H 2 ) separation 



C0 2 -permeability 

9.0 (5.6) barrer 


Selectivity C0 2 /CH 4 , C0 2 /C0, C0 2 /N 2 

21.1 (22.4) 


Selectivity H 2 /CH 4 

6.2 (56.0) 


Effective membrane thickness 

1000 A 

All sections 

Heat exchanger pressure drop 

neglected 


Isentropic efficiency of turbomachinery 

80% 

Thermodynamic models 

Ideal gas law, liquid phase in unsymmetric convention with solubilities from [52] 
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a more conservative hypothesis concerning methane forma¬ 
tion in the model, which is to consider not the partial pres¬ 
sure, but the molar fractions as activity of the species in the 
gas phase. A K~ is therefore defined, i.e.: 




K e=ff J* = K fP j 


-(E"rE n i> 


(17) 


ner 

i 

and used instead of K p . Equation (14) is therefore rewritten as: 


K e = K p (T g +AT eq ) (18) 

in which K p is still evaluated with Equation (15). 

By applying this approach to four independent equilibrium 
reactions of Table 2, it is feasible to determine the reaction 
extent and the product composition of gasification. However, 
it has been observed that the observed C 2 H 4 -fraction is 
excessively far from its equilibrium composition in terms of 
temperature, and that the amount of residual carbon is 
numerically sensible and difficult to handle since it is only 
implicitly represented in the equations. In order to avoid 
a considerable deterioration of the robustness of the model, 
the corrected equilibrium reactions are therefore only used 
with the hydrogenating gasification (3) and the water-gas shift 
reaction (7). The amount of higher hydrocarbons represented 
by ethene is assumed to be proportional to the methane yield 
(19), and the carbon conversion efficiency is considered 
constant (20): 


Cc 2 H 4 — kc 2 /CiCcH 4 


(19) 


^■carbon,residual — (1 ^cc) ^-carbon, biogenic 


( 20 ) 


Table 5 - 

Reconciled gasification model parameters. 

Gasifier 

Indirectly heated 

Directly heated 

AT eq,(3) 

-280 °C 

-280 °C 

AT eq) ( 7 ) 

-112°C 

-112°C 

k C 2 /Ci 

0.205 

0.476 

£ C c 

90.3% 

93.0% 


of the model parameters of Table 5. It is thereby interesting 
to see that the distance to equilibrium of the reactions are 
identical in the two models although the principles of gasi¬ 
fication are different. This means that the bed plays the same 
catalytical role in the gasifier and that the stability of 
methane in the gas phase preventing its decomposition is the 
same in the two gasifiers which are operating at similar 
temperatures. 

Fig. 3 depicts the dependence of the gasifier outlet 
composition in the interval of ±50 °C of the nominal gasifi¬ 
cation temperature. For both gasifiers, lower operating 
temperatures are increasing the hydrocarbon and C0 2 content 
and decreasing the H 2 and CO fractions, resulting in an 
increased volumetric calorific value of the gas. Comparing 
both gasifiers, it can be seen that the CO and C0 2 fractions and 
thus the degree of oxidation is considerably higher in directly 
heated gasification. Accordingly, its product gas has a lower 
stoichiometric number than in indirectly heated gasification, 
which will lead to a lower CH 4 concentration after methana- 
tion. Fig. 4 shows the cold gas efficiency £ cg , defined as the ratio 
of the chemical energy contained in the cold product gas and 
the raw material (21), for variations of temperature and inlet 
wood humidity around the normal operating point. 


By fitting the model parameters (AT eq> ( 2 ), AT eq> ( 2 ), fe^/ca and 
£ cc ) to data of existing plants, both correct product yields and 
energy balances around the nominal operating points are 
obtained. Table 4 shows the good agreement of the recon¬ 
ciliation reached for the considered gasifiers with the values 


Ah Lod™ 


gas 

+ 

wood 


( 21 ) 


According to the figure, the conversion of chemical energy is 
more efficient in directly heated gasification and situated 
around 80%, while its value for indirectly heated gasification is 
in the range of 70-80%. As discussed in Section 2.2.1, 


Table 4 - Reconciliation of the producer gas compositions 


at nominal T g and p g = 1 bar (Data/Calculation [%vol]). 


Gasifier 

Indirectly 
heated [66] 

Directly 
heated [29] 

Directly 
heated (0 2 ) 


850 °C 

800 °C 

800 °C 

C 2 H 4 

1.8/1.9 

2.0/1.9 

74.5 

ch 4 

8.8/9.6 

4.2/4.0 

79.5 

h 2 

37.3/38.5 

14.8/14.7 

725.8 

CO 

29.4/27.4 

15.4/16.0 

732.3 

co 2 

16.2/15.8 

15.0/14.7 

724.0 

n 2 

-72.9 

39.6/40.3 

70.1 

h 2 o 

3.6/3.9 

78.4 

73.8 

Ah 0 

12.0/12.2 

6.2/6.2 

712.9 

[MJ-Nm -3 ] 
SN [-] 

0.24/0.26 

0.13/0.13 

70.13 


a Although no nitrogen is introduced by the gasification agent, 
some N 2 is used for inertisation of the raw material, which 
prohibits to attain the criterion on the Wobbe Index at the process 
outlet. In the remainder of this work, a cut-down to 0.5%vol of the 
dry gas by inertisation with C0 2 is assumed feasible. 



Gasification temperature [°C] 

Fig. 3 - Molar compositions of the dry gas from gasification. 
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Fig. 4 - Cold gas efficiencies in parameter domain. 


increasing wood humidity is markedly deteriorating the 
performance of both gasifiers. 

3.1.3. Oxygen production for gasification 

Since no major advantages from the integration of the SNG 
production with a cryogenic or adsorptive air separation plant 
is expected, oxygen has been considered as a utility and the 
cost figures from Kirschner [30] are used directly. If by- 
produced oxygen from electrolysis is used, an energy effi¬ 
ciency of 85% based on the lower heating value of hydrogen 
and the power input is assumed, and the excess heat is 
considered available at 120 °C. A detailed analysis of this 
particular process configuration is reported in [32] and thus 
not discussed in this paper. 

3.1.4. Thermal pretreatment before gasification 
Compared to gasification, torrefaction and pyrolysis are 
processes where the solid decomposition is less advanced and 
multiple condensable and non-condensable products are 
formed. An equilibrium approach to predict accurate product 
yields and composition is not adequate, and the processes are 
represented by simple conversion ratios for fixed operating 
conditions. From experimental torrefaction data at 260 °C [67], 
a dry solid yield of 87% has been derived, whereas the carbon 
mass fraction in the product is increased by 5.7% compared to 
its initial hydrocarbon composition, and the hydrogen and 
oxygen fractions are decreased by 8.4% and 5.8%, respectively. 
The moisture in the vapour phase is accompanied by acetic 
acid, methane, carbon mono- and dioxide, whereas 50.1% and 
17.5% of the biogenic hydrogen and oxygen is bound in 
methane and carbon monoxide, respectively. For screw 
pyrolysis [36], it is assumed that the volatile fraction of the 
solid is completely transferred into the vapour phase, and 
a simplified model for calculating the heat requirement for the 
decomposition is used [22]. 

3.1.5. Gas cleaning 

Since the conversion of impurities and tars are negligible for 
energy concerns, they have not been included in the models. 


However, the thermodynamic transformations that the 
cleaning operations imply obviously need to be considered. If 
cold cleaning is applied, the producer gas from gasification is 
cooled to 150 °C before entering the filter. Downstream of this 
unit, no thermal energy can be recovered and the gas is cooled 
to atmospheric temperature in the scrubber, where it also 
reaches water saturation. In the model, it has been assumed 
that this heat is lost, although part of it would be available as 
moderately heated water leaving the scrubber. If hot cleaning 
is applied, the sensible heat of the gas is supposed to be 
entirely recoverable. 

3.1.6. Methane synthesis 

Using a similar modelling approach like in gasification, the 
outlet composition of the methane synthesis reactor has 
been computed with equilibrium equations. However, 
methanation is a catalytic process and the equilibrium 
condition is a reasonable assumption if the amount of cata¬ 
lytic material is sufficient. Experimental data from a labora¬ 
tory reactor confirm this assumption [5]. No model constants 
have therefore been introduced and the methanation (5), the 
steam ethene reforming (6) and the water-gas shift reaction 
(7) have been supposed to be in chemical equilibrium. In 
order to avoid carbon deposition, steam must be added to the 
producer gas prior to the reactor which disfavours the 
formation of methane. The minimum required amount is 
estimated based on the equilibrium data presented by 
Mozaffarian and Zwart [3]. 

3.1.7. Carbon dioxide removal 

All considered carbon dioxide removal processes are based on 
local diffusion processes and are heavily dependent on 
diffusion and absorption kinetics. A detailed model requires 
dynamic simulation, which is too complicated and not 
appropriate for flowsheet calculations in process design 
studies. For this reason, the C0 2 -removal models are devel¬ 
oped on the basis of overall performance and characteristic 
operation parameters. 

For physical absorption, the solubilities of carbon dioxide 
and methane in the classic solvent Selexol reported by Sweny 
and Valentine [48] are used to calculate the multicomponent 
separation in an absorption tower. Applying the Kremser 
method [68], its performance is determined with respect to the 
number of theoretical trays and the relative solvent flow rate. 
While the residual humidity in the saturated gas is completely 
removed by the solvent, the solubility of the other non¬ 
condensable species (H 2 , CO, N 2 ) is very low and therefore 
neglected in the modelling. In order to achieve a high product 
recovery, a typical process layout with a flash drum at 
moderate pressure for recycling the dissolved CH 4 is consid¬ 
ered [69]. 

Pressure swing adsorption is a discontinuous process that 
removes the carbon dioxide by its adsorption under pressure 
following regeneration of the adsorbent at subatmospheric 
pressure. The purity and the amount of methane recovered 
in the outlet stream is essentially determined by the dura¬ 
tions of the adsorption, recycling and purging periods. As 
shown on Fig. 5, two parameters, i.e. t n and t r2 , are intro¬ 
duced and fix the relative durations of these periods. The 
time-averaged flow of species i that leaves the adsorber 
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C0 4 0.0058 0.0305 -0.2957 1.2415 -0.6652 


CH 4 and others 0.8148 0.3758 0.2021 -0.8486 0.4547 

Fig. 5 - Illustration of the phenomenological model for PSA 
and coefficients for the regression of f(i, t r ) = a 0 + a t t r 
+ a 2 t? + a 3 t? + c* 4 t y used in Equations (22) and (23). 


system (m iout ) or is recycled to its inlet (rh irec ) is then deter¬ 
mined by a relation of the form: 


m im =/(i,t r i)m y „ 

(22) 

iili.rec =f(ht' r2 )m iM -m iout 
with t' r2 = (1 - tn)t r2 + tn 

(23) 


where / (i, t r ) represents the relative amount of species i that 
has left an adsorption vessel at t r during one cycle. A numeric 
integration of pilot plant data reported by [49] for biogas 
separation with a carbon molecular sieve adsorbent has 
allowed for regressing / (i, t r ) to a polynomial of 4th degree 
whose coefficients are shown on Fig. 5. 

For SNG-upgrade with membranes, a general design model 
for hollow-fibre modules proposed by [70] has been imple¬ 
mented in a superstructure of possible membrane arrange¬ 
ments. At the example of this separation technology, the 
integration of the crude SNG separation with the reactive 
process steps is discussed and optimised in a separate work 
[71]. As also for the PSA process, it is assumed that the satu¬ 
rated feed is first dried by temperature swing adsorption (TSA) 
over aluminium-oxide, for which HMJ-kg^o at 160-190 °C 
are required [72]. 

Exploring these models, it has been determined that the 
Wobbe Index of 13.3 kWh-Nm 3 can be met with all considered 
technologies. For upgrading the crude gas from atmospheric 
pressure to grid quality and pressure by PSA or physical 
absorption, a specific mechanical power consumption in the 
order of 600-800 kj kg^i deS NG is required and leads to an 
approximate composition of 84%vol CH 4 ,12% H 2 , 0.2% CO and 
3.5% of inerts. With these two technologies, the residual 
hydrogen from the equilibrium synthesis is however not 
separated from methane, which limits its purity to around 
86%vol CH 4 . For gas upgrading to a methane content of 
96%vol, a final membrane stage for hydrogen removal can be 
used, from which the H 2 -rich permeate is recycled to the 
methane synthesis. According to the information provided by 
Abetz et al. [73], matrimid or polysulfone membranes are 
a reasonable choice for this purpose. Alternatively, an 
upgrading system entirely based on cellulose acetate 
membranes can be used, for which a specific mechanical 
energy consumption of 750 to 1100 kj kg^ r i deSNG depending on 
the targeted SNG-recovery in the separation stage has been 


assessed [71]. In case of indirectly heated gasification, the 
depleted stream from the gas upgrading is thereby efficiently 
used as fuel, which allows for limiting the size of the separa¬ 
tion system due to a better process integration. 

3.2. Energy-integration model 

As discussed in detail in the applied methodology [15], the 
heat requirements of all process stream determined in the 
energy-flow models are implemented in the energy-inte¬ 
gration model to compute the optimal thermal integration of 
the process streams and to define the heat recovery poten¬ 
tial. First, the temperatures of all hot and cold streams are 
corrected by minimum approach temperatures that depend 
on the type of stream considered 2 . Then, the temperature 
intervals assembling the streams are identified and the 
energy cascade of the process is calculated. In general, 
several possibilities to satisfy the minimum energy 
requirements are provided. The necessary heat that needs to 
be supplied above the pinch is usually produced by 
combustion. If available at a useful temperature, the excess 
heat below the pinch is recovered in a Rankine cycle for 
power production, whose superstructure (i.e. the production 
and usage levels) is predefined in the problem set-up. In this 
work, a steam cycle with back pressure and condensing 
turbine stages at efficiencies of 80% and 70%, respectively, 
are considered. Depending on the possibility of selling heat 
in a district heating network, the cycle may be designed for 
power- or heat-operation. Other utility techniques, like the 
use of heat pumps, combustion boosting through oxygen 
injection or air preheating might be used if large quantities 
of heat are exchanged around the pinch point or if the pinch 
is situated at a very high temperature level. 

A mathematical programming formulation of the optimal 
integration of utilities has been proposed and detailed in 
[17,15]. The method allows for computing the optimal flows 
through the system for maximum combined fuel, power and 
heat production. In this approach, we consider off-gases and 
depleted streams from the different process sections as fuel 
and balance the requirement by withdrawing intermediate 
product streams from suitable locations. The flow rates of the 
different sections are thus interdependent with the heating 
and cooling requirements of the process, and the optimal 
system is determined by mixed integer linear programming 
maximising the operating profit from the products. 

3.3. Economic evaluation 

3.3.1. Approach 

In the design methodology applied here, we consider the 
thermodynamic conditions to be reached as decision vari¬ 
ables. Once these conditions are determined in the flowsheet 
model, we compute the process performances and estimate 
its investment cost by sizing the major process equipment 
that allows to reach the thermodynamic design target. This 
method differs from the conventional thermo-economic 


2 Minimum approach temperature contributions of 8, 4 and 2 K 
for gaseous, liquid and evaporating or condensing streams are 

assumed here. 
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approaches where either the equipment size is considered as 
a decision variable or the total investment cost is estimated on 
the components’ nominal capacity, as it is done by for 
example [8]. 

The investment estimation presented here is obtained by 
realizing a simplified sizing of the major process equipment. 
Considering that the equipment’s cost is essentially dependent 
on its size and construction materials, it is influenced by 
specific process conditions like temperature, pressure and 
volume flows. Hence, the information available from flowsheet 
calculations allow to relate the process investment cost with its 
thermodynamic operating conditions and performances. This 
is of great interest with regard to a future thermo-economic 
optimisation of the installation. Following the method 
described in [19], the total grass roots costs C GR that designate 
the total investment cost for a new production site, can be 
correlated to the equipment’s purchase cost by (24): 


Cgr - (1 + Cl) ^ C BM ,i + c 2 ^ j (24) 

i i 

where C° BM) i is the bare module costs of the equipment i at 
base case conditions (i.e. ambient operating pressure and 
carbon steel construction) and C B M,i its cost considering the 
actual operating conditions. The factors Ci and c 2 represent 
additional expenses related to the construction of the plant. 
Following Turton et al. [19], it is assumed that the costs of 
auxiliary facilities, site development and buildings are inde¬ 
pendent on special process conditions and amount to 35% of 
the total bare module costs at base case conditions. Contin¬ 
gencies and fees are dependent on these conditions and 
amount to 18% of the equipments’ total cost. The according 
numeric values of 0.18 for c 1 and 0.35 for c 2 used here provide 
an investment cost estimation for mature technology and 
established process engineering. 

Unless stated otherwise, the bare module costs of the 
equipment is determined from correlations given in [18,19]. 
For process vessels, maximum diameters of 4 m (vertical) and 
3 m (horizontal) are assumed, and multiple units operated in 
parallel are considered if necessary. A complete list of 
parameters for the investment analysis and the considered 
prices for wood and electricity are given in Table 6. 


3.3.2. Equipment rating 


3.3.2.1. Wood Drying. The sizing of the direct rotary air dryer 
is based on an average velocity of wood of u wood = 0.03 m s -1 
occupying 12% of the dryer’s cross sectional area [18]. 
Accordingly, the diameter of the drum is given by: 


d = 2 


m wood 


0.12I7p wood u wood 


1/2 


(25) 


where rh wood and p WOO d are the mass flow and density of wood 
respectively. The dryer’s length is determined considering an 
overall heat transfer coefficient in W m _3 K _1 based on the 
dryers volume as given by [18]: 


U = 240G 067 /d (26) 

where G is the average gas mass flux that is typically in the 
range of 0.5 to 5 kgs -1 m -2 . The dryer’s length follows from: 


Ah va pVh]-i 2 

O ,vap ( r )7\ 

UAAT lm { } 

For steam drying fixed bed design is a convenient choice and 
detailed information about its performance and size is given 
in [25]. Assuming the mean velocity of superheated steam 
being constant during scale-up, the dryer’s diameter is 
calculated with: 


d = 2 


( V s team 
nUynean 


1/2 


(28) 


where V steam is the volume flow of steam and u mean its mean 
velocity. Using the developed energy-flow model for steam 
drying and data from [25], a value of 1.4ms -1 has been 
assessed for the latter. The height of the steam dryer is 
determined in the same way as for air drying and Equations 
(26) and (27) are used to calculate the overall heat transfer 
coefficient and the height of the bed. Thereby, the correlation 
(Eq. (26)) has been reconciled with the data given in [25] and 
very good accordance is observed. At an average gas mass flux 
of 0.645 kg s -1 m -2 , the observed value for U amounts to 
200 W m _3 K _1 , while its calculated value is of 199 W m -3 K -1 . 
Finally, taking the additional height of the dryer’s freeboard 
above the bed into account and preserving the geometry of the 
considered steam dryer, the overall height of the vessel has 
been obtained by multiplying the bed height with a factor of 
1.27 [25]. For both the direct rotary air dryer and the fixed bed 
steam dryer, costing data from [18] have been used. 


3.3.2.2. Gasification and methanation. The gasification and 
methanation reactors are of fluidised bed type and, as it is the 
case for the steam dryer, their mean gas velocities will typi¬ 
cally remain constant during scale-up. The diameter of the 
reactors are thus calculated according to Equation (28), 
replacing V s team by the total gas volume flow V gas . To calculate 
the reactor heights, an exponential law of the form: 

h = hoC (29) 


Table 6 - Assumptions for process 

economics. 

Parameter 

Value 

Marshall and Swift index (2006) 

1302 

Dollar exchange rate 

1 €/US$ 

Interest rate 

6% 

Discount period 

15 years 

Plant availability 

90% 

Operators 3 

4 b p./shift 

Operator salary 

60’000€/year 

Maintenance costs 

5%/year of C GR 

Wood price (<t>u; 00d = 50%) 

33 €/MWh 

Biodiesel price 

105 €/MWh 

Electricity price (green) 

180 €/MWh 


a Full time operation requires three shifts per day. With a working 
time of five days per week and 48 weeks per year, one operator per 
shift corresponds to 4.56 employees. 

b For a plant size of 20 MW th|UJOOd . For other production scales, an 
exponent of 0.7 with respect to plant capacity is used. 
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has been fitted to data from existing plants [46,74]. The 
correlated data for Equations (28) and (29) are given in Table 7. 

In order to take the costs associated with the gasification 
reactor internals and their special construction into account, 
the costs of the vessels have been multiplied by a factor of 4.4. 
This value has been derived from the total grass roots costs of 
the existing FICFB pilot plant, which is of 8 M€ for an instal¬ 
lation with a nominal capacity of 8MW th [75]. In case of 
a pressurised gasification, the investment for the biomass 
feeding system may become substantial and is thus added 
explicitly. The use of lock hoppers is the most common 
practice [76], and a manufacturer’s cost estimate reported by 
Swanson et al. [[Appendix G] [77]] is used to determine 
a standard scaling law for a double-train system: 

Cbm, lock hoppers =4.3-10 -VYl drywood (30) 

with Cbm, lock hoppers in US$ (2002). 

For oxygen supply to directly heated gasification, cost data 
from [30] are used, who reports a range of 0.03-0.70 US$kg -1 
(1999) depending on the quantity required. In case production 
with electrolysis is considered, the bare module cost of the 
unit is assumed according to the US/DOE target (2010) of 
300 $ kWej 1 [78]. 

3.3.2.3. Thermal pretreatment before gasification. According 
to a commercial design [79], the torrefaction reactor is 
dimensioned as a vertical tower with a height to diameter 
ratio of 2, assuming that 20% of the volume are occupied by 
the solid with a residence time of 15 min. The purchase cost of 
the unit has been regressed on the reported 5.65-10 5 and 
6.70-10 5 US$ for units with a volume of 115 and 145 m 3 , 
respectively. 

For pyrolysis, the design of Henriksen et al. [36] has been 
investigated and procedures for a screw conveyor or rotary 
calcinator are used [18], where the unit size is determined by 
the required heat transfer through the wall to the solid. With 
a heat transfer coefficient of 30-100 kW m -2 K -1 , scale-up of 
the pilot plant is however expected to be costly, since either 
a large temperature difference or heat transfer area is needed. 

3.3.2.4. Gas cleaning. While the cost of a cold gas cleaning 
system including a cyclone, bag filter, scrubber and guard 
beds can be estimated from literature data [18], hot gas 
cleaning technology is still under development and no cost 
data is currently available. It is therefore assumed that 
mature hot gas cleaning will roughly follow the same 
correlation than cold technology, and its cost is evaluated 
with the same relations using the gas volume flow as 
parameter. At an assumed temperature of 500 °C, the gas 


volume flow to be treated is nearly doubled compared to 
cold processing, which leads to a multiplication of the 
investment in the order of 1.5. 

3.3.2.5. Carbon dioxide removal. The dimensions of the 
packed tower in the physical absorption process have been 
calculated according to the method outlined in [18]. From the 
given data, a stage efficiency of 15% has been determined, and 
the tower diameter has been dimensioned considering the 
entrainment limit. The corresponding stripper has not been 
modelled in detail. As a rule of thumb, its diameter and height 
have been assumed to be equal to the ones of the absorption 
tower [48]. 

The size of the four pressure swing adsorption vessels 
necessary for continuous operation has been estimated 
assuming that the required adsorber volume is proportional to 
the amount of carbon dioxide that needs to be removed from 
the gas. A typical plant treating 1000 Nm 3 raw gas per hour 
with a C0 2 content of 40%vol is described in [80]. The diameter 
and height of the vessels are reported to be 1.4 and 6 meters, 
resulting in an specific adsorber volume of 83.13 m 3 per 
Nm^s -1 . The shape of the vessels has been assumed 
constant and a height to diameter ratio of 4.3 is used. 

The purchase cost of the polymeric hollow-fibre 
membranes is estimated by updating the data reported by 
Bhide and Stern [50], who considered an initial cost for the 
permeator modules of 108 US$ m -2 and a membrane element 
cost of 54 US$ m -2 that need to be replaced every three years. 
As detailed in [71], the required membrane area for the 
separation is directly computed by the design model. 

In case of separation by PSA or membranes, the cost for the 
required TSA unit used for gas drying is determined by 
assuming a two-column layout with a cycle time of 12 h and 
a maximum adsorbent loading of 0.12 kg^-kg^,.^ [81]. 
The cost of the adsorbent and its density have been assumed 
to 9 US$kg -1 and 800 kgm -3 , respectively [18]. 

If C0 2 -removal is done prior to methanation, the cost esti¬ 
mation of the necessary shift reactor to adjust the stoichiom¬ 
etry is based on the preliminary design and data from [13]. 

3.3.2.6. Heat exchangers and turbomachinery. The total cost 
of the heat recovery system including the heat exchangers of 
the process and the utility system has been computed from 
the hot and cold composite curves assuming an overall heat 
transfer coefficient of 580 W m -2 K -1 . The total heat exchange 
area is distributed over the minimum number of heat 
exchangers and costing data for fixed tube sheet units from 
[18] have been used. For turbomachinery, data for mainly 
centrifugal units are used. 


Table 7 - Sizing parameters for process reactors. 
b = 0.188 is reconciled for the gasification vessel and 
used for all reactors. Data calculated from [46,74]. 

Reactor type 

Mmean [^n S ] 

h 0 [m-(m 3 -s *) b ] 

Gasification 

0.645 

4.07 

FICFB-combustion 

5.250 

8.47 

Methanation 

0.093 

18.0 


4. Performance of exemplary technological 
options 

The outlined thermo-economic model is developed for 
a systematic process optimisation, which exceeds the scope of 
this paper and will be presented separately. Here, the 
modelling approach is thus only illustrated at the example of 
a few technology scenarios shown in Table 8, which should 
not be considered as a portfolio of optimal process options. 
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Indirectly heated gasification with air drying and PSA for 
removing C0 2 after methanation is defined as a base scenario 
(FICFB, base), to which different alternatives are compared. In 
the scenario (torr), a torrefaction unit is added to the base 
case, (pM) investigates the influence of a moderately pres¬ 
surised methanation, and (pM, SA) compares the latter case to 
one where the stoichiometric number is adjusted to 1 by 
water-gas shift conversion and C0 2 -removal prior to metha¬ 
nation. Indirectly heated gasification is then compared to 
pressurised steam/oxygen-blown gasification and methana¬ 
tion at the same pressure (CFB, pGM), for which the potential 
improvement by introducing hot gas cleaning is assessed 
(pGM, hot). In all cases, excess heat is recovered in a steam 
Rankine cycle, whose production and utilisation levels have 
been adjusted to reasonably fit the steam demand for gasifi¬ 
cation and methanation, and providing heat for mainly drying 
and torrefaction. Heat extraction for district heating is not 
considered, and a condensation level suitable for the cold 
utility is used. 

4.1. Energetic performance 

Table 9 shows the overall energy balance and conversion 
efficiencies of the scenarios defined in Table 8 for a capacity 
of 20 MW th>wood . One integrated composite curves for each 
gasification technology is given in Fig. 6, and Table 10 shows 
the stream compositions and flows through the systems. 
Although the energy balances of the different cases seem 
rather close, some general trends depending on the tech¬ 
nology choice can be observed. Due to the higher cold gas 
efficiency for directly heated gasification assessed in Fig. 4, 
the resulting SNG-yield is roughly 10% higher than in case of 


indirectly heated gasification. As shown on the composite 
curve, the indirect heat supply creates a process pinch point 
at the gasification temperature, and satisfying this heat 
demand at high temperature requires to withdraw about 
19% of the raw gas production of the gasifier. For this 
reason, more excess heat is also available below the pinch 
and leads to a higher co-production of electricity than in 
directly heated gasification. An alternative to using the 
excess heat in a Rankine cycle is thereby to recover it in the 
process through thermal pretreatment. Acting in principle 
as a chemical heat pump, the endothermal decomposition 
starts with excess heat from below the pinch, and less heat 
must be supplied above the process pinch point [22]. In case 
of torrefaction, it is necessary to withdraw only 5% of the 
raw gas production for satisfying the process energy 
requirement, which increases the overall SNG-yield from 
the indirectly heated gasifier by 7%. 

While the methanation pressure and prior stoichiometry 
adjustment influences the energetic performance to a lesser 
extent, hot gas cleaning is expected to considerable increase 
the process performance by cogenerating more electricity. 
Heat losses due to cooling, washing and reheating the gas are 
avoided, and less steam needs to be prepared for the metha¬ 
nation. More heat can thus be recovered in the Rankine cycle, 
whose gross production increase by 0.5 MW e j and results in an 
positive electricity balance of the overall system also in case of 
directly heated gasification with oxygen. 

4.2. Economic performance 

With a common share of 68% of the total investment of 32.6 
M€ illustrated in Fig. 7, the gasification and methanation 
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Table 9 - Useful energy balance and efficiencies of the overall system. 


FICFB CFB 





(base) 

(torr) 

(pM) 

(pM, SA) 

(pGM) 

(pGM, hot) 

Consumption 

Wood 

kW 

20’000 

20’000 

20’000 

20’000 

20’000 

20’000 


Biodiesel 

kW 

365 

316 

365 

365 

22 

- 


Electricity 

kW 

- 

63 

- 

- 

46 a 

- 

Production 

SNG 

kW 

13’443 

14’318 

13’446 

13’495 

14765 

14769 


Electricity 

kW 

609 

- 

521 

526 

- 

315 a 

Efficiency 13 

£ 

% 

69.0 

70.3 

68.6 

68.8 

73.6 

75.4 


c 

^chem 

% 

71.3 

70.1 

70.5 

70.8 

73.5 

76.6 


V 

% 

63.5 

64.7 

63.1 

63.3 

67.8 

69.4 


a A consumption of 311 kW (1080 kj-kg^) for off-site oxygen production is included, 
b Energy e and exergy rj efficiencies are defined as the ratio between production and consumption terms. 

c Chemical efficiency: electricity is expressed as SNG-equivalent and substituted therein. An exergy efficiency of 55% is assumed for the 
conversion of SNG to electricity. 


sections are dominating the capital costs of the (base)- 
scenario for FICFB-gasification. As shown in Table 11, atmo¬ 
spheric pressure leads to big reactors due to relatively 
important volumetric flows, and a much more economic plant 
design can be obtained by pressurising the reactors. Already at 
a moderate methanation pressure of 5.5 bar, the investment 
cost are predicted to decrease by 30% to 23.4 M€. If directly 
heated gasification with oxygen is used and the whole system 
is initially pressurised to 15 bar, a further reduction to 17.9 M€ 
is possible, whereas the pressurised feed train accounts for 
approximately 10% of the gasifier cost. If a scale-up to larger 
plant capacities is considered, operation at atmospheric 
pressure requires several units in parallel, and a pressurised 
system seems unavoidable. 

Based on the assumptions of Table 6, the share of the 
depreciated investment and its related maintenance cost is 
relatively high. At base case conditions, these expenses 
amount to 47% of the total production cost. If methanation 
is moderately pressurised, the expected reduction of the 
investment allows for decreasing the share of these 
expenses to 38% of the total, and leads to a cost reduction by 
13% in absolute terms. In the best performing case of 



directly heated gasification with hot gas cleaning, the share 
of investment and maintenance is of only 32% at 
20 MW th)U ;ood, and can further be reduced to 23% in a plant 
scale-up to 150 MW th>wood . For this directly heated option, 
the cost for oxygen is thereby not excessively penalising the 
plant economics. With the cost data from Kirschner [30], the 
required 0.29kgs -1 (820Nm 3 h _1 ) at 20MW th)U;ood are best 
supplied by on-site PSA. At this scale, a price of 9.4 cts kg^ 
is obtained, and the expenses for oxygen amount to 
4.4 € MWh _1 SNG , which correspond to 5-6% of the total 
production costs. 

Overall, directly heated gasification tends to be the better 
option, for which total production costs of 80.4 
(75.7) €MWh -1 S NG with cold (hot) gas cleaning are assessed. 
For a plant at 150 MW th)WOod , the total costs would reduce to 
63.6 (58.9) € MWh -1 SNG . These numbers are by 10-20% better 
than in indirectly heated gasification with moderately pres¬ 
surised methanation, for which the production costs amount 
to 90.9 and 80.1 € MWh -1 SN G at 20 and 150 MW thtWOod , respec¬ 
tively. Compared with an expected market price of 
120 € MWh -1 S NG for the produced green gas, all these options 
are yet expected to be profitable. 



Fig. 6 - Integrated composite curves for FICFB (base, on the left) and CFB (pGM,hot) for an input of 20 MW th>WOO d. 
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Table 10 - Dry gas composition, calorific value and power balance at process section outlets. 


Gasification 

Methanation 

SNG-upgrade 

FICFB CFB 

FICFB a CFB 

FICFB a CFB 


c 2 h 4 

%vol 

2.3 

5.5 

- 


- 

- 

(-) 

- 

ch 4 

%vol 

11.3 

11.6 

48.3 

(86.1) 

40.8 

96.0 

(96.0) 

96.0 

h 2 

%vol 

40.2 

28.9 

6.4 

(11.4) 

1.4 

1.6 

(1.9) 

1.3 

CO 

%vol 

26.5 

17.4 

0.1 

(0.0) 

0.0 

0.3 

(0.0) 

0.1 

co 2 

%vol 

19.2 

36.5 

44.4 

(1.1) 

57.6 

0.5 

(0.5) 

2.2 

n 2 

%vol 

0.5 

0.1 

0.8 

(1.4) 

0.2 

1.6 

(1.6) 

0.4 

Ah 0 

MW kg -1 

15.2 

11.4 

14.5 

(47.4) 

10.3 

48.2 

(47.6) 

46.8 

Wobbe Index 

kWh Nm“ 3 

4.3 

4.3 

5.6 

(13.6) 

4.3 

14.1 

(14.1) 

13.9 

Load 

MW th 

19.6 

17.0 

15.1 

(14.7) 

15.6 

13.4 

(13.5) 

14.8 


a Values in parenthesis are for C0 2 -removal before methanation in case (pM, SA). 
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FICFB 


(pM,SA) ; (pGM) (pGM,hot) 
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Fig. 7 - Investment (left) and total production costs for the scenarios defined in Table 8 at a plant capacity of 20 MW t h )U , 0 od« 
Negative contributions are due to profits from selling electricity. 


Table 11 - Reactor size and number of required units for 
a plant capacity of 20 MW thiU , 0 od (150 MW th>luood , 
dimensions change). 


FICFB 


CFB 


(base) (pM) 

(pM,SA) 

(pGM) 

N 3 

1 (8) 1 (8) 

1(8) 

1(1) 

dg/ ^g,comb. 

3.9 m/1.1 m 


1.0 m 

hg/hg,comb. 

6.0 m/11.6 m 


3.5 m 

N m 

3 (21) 1 (4) 

1(3) 

1(2) 

d m 

3.8 m 2.8 m 

2.5 m 

1.9 m 

K 

18.2 m 16.2 m 

15.4 m 

14.0 m 


5. Conclusions 

This paper presented a superstructure based process model 
for candidate technologies to produce SNG from woody 
biomass. Based on data from existing plants and pilot instal¬ 
lations, simple parametric flowsheet models have been vali¬ 
dated and are expected to accurately predict the 
performances of the installation in a restricted domain of 
operating conditions. Compared to conventional flowsheet 


simulation, the proposed model has been developed to assess 
the thermo-economic performances of the integrated system. 
It thereby computes not only the energy-flow using conven¬ 
tional flowsheeting software, but also the process integration 
including the combined heat and power production. The 
energy-integration is performed without imposing any 
restrictions on the process topology, and therefore expected to 
determine the optimal choice of utilities. 

The results of this stepwise flowsheeting are then used to 
rate the equipment and estimate the investment cost of the 
system. For these economic calculations, thermodynamic 
variables like temperature, pressure and volume flows are 
considered as design specifications to be met, which allows to 
take the impact of the operating conditions into account. Due 
to this structure, the process model is adequate for a future 
thermo-economic optimisation. 

In a preliminary screening of the process performance 
for some typical technology scenarios, it is shown that the 
conversion of woody biomass to SNG is a viable option with 
respect to both energetic and economic aspects. Recovering 
the excess heat by means of a Rankine cycle allows for 
a considerable co-production of electricity, and overall 
energy and exergy efficiencies of in the range of 69-76% and 
63-70%, respectively, are obtained. The most efficient 
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conversion is thereby reached with directly heated, steam/ 
oxygen-blown gasification and hot gas cleaning. Due to 
reduced volume flows, moderately pressurised gasification 
and methanation is further advantageous with regard to 
investment cost. Including the equipment’s depreciation, 
total production costs in the range of 76-107 € MWh _1 SN G 
are expected for a plant capacity of 20 MW th>wood , whereas 
59-97 € MWh _1 SNG are assessed for large-scale plants at 
150 MW thiUJOO£ j and above. 
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